Draft version February 1, 2008 

Preprint typeset using LAT^^X style cmulateapj v. 14/09/00 



RADIO ASTRONOMICAL POLARIMETRY AND PHASE- COHERENT MATRIX CONVOLUTION 

W. VAN StRATEN 

Centre for Astrophysics and Supercomputing, Swinburne University of Technology, 
Hawthorn, VIC 3122, AustraUa 
wvanstraCipulsar . physics . swin.edu.au 
Draft version February 1, 2008 

ABSTRACT 

A new phase-coherent technique for the cahbration of polarimetric data is presented. Similar to the 
one-dimensional form of convolution, data are multiplied by the response function in the frequency 
domain. Therefore, the system response may be corrected with arbitrarily high spectral resolution, 
effectively treating the problem of bandwidth depolarization. As well, the original temporal resolution 
of the data is retained. The method is therefore particularly useful in the study of radio pulsars, where 
high time resolution and polarization purity are essential requirements of high-precision timing. As 
a demonstration of the technique, it is applied to full-polarization baseband recordings of the nearby 
millisecond pulsar, PSR J0437-4715. 

Subject headings: methods: data analysis — pulsars: individual (PSR J0437-4715) — techniques: 
polarimetric 



1. INTRODUCTION 

In radio polarimetry, two orthogonal senses of polariza- 
tion are received and propagated through separate signal 
paths. Each signal therefore experiences a different series 
of amplification, attenuation, mixing, and filtering before 
sampling or detection is performed. Whereas efforts are 
made to match the components of the observatory equip- 
ment, each will realistically have a unique frequency re- 
sponse to the input signal. Even a simple mismatch in 
signal path length will result in a relative phase difference 
between the two polarizations that varies linearly with fre- 
quency. 

In fact, any physically realizable system will transform 
the radiation in a frequency-dependent manner. Where 
variations across the smallest bandwidth available may be 
considered negligible, post-detection calibration and cor- 
rection techniques may be used to invert the transforma- 
tion and recover the original polarimetric state. However, 
the transformation may vary significantly across the band, 
causing the polarization vector to combine destructively 
when integrated in frequency. This phenomenon is known 
as "bandwidth depolarization" of the signal, and results 
in irreversible decimation of the degree of polarization. 

It is therefore desirable to perform polarimetric correc- 
tions at sufficiently high spectral resolution, which is avail- 
able only at the cost of temporal resolution unless phase 
coherence is maintained. Conventional post-detection cor- 
rection techniques therefore prove insufficient when high 
time resolution is also a necessity. For example, certain 
pulsar experiments require high time resolution in order to 
resolve key features in the average pulse profile. As well, it 
has been shown that insufficient time resolution can also 
lead to depolarization (Gangadhara et al. 1999). These 
considerations motivate the development of a method for 
phase-coherent polarimetric transformation. Further im- 
petus is provided by the growing number of baseband 
recording systems at observatories around the world, the 
enhanced flexibility made available through use of off-line 



data reduction software, and the increasing computational 
power of affordable facilities. 

A method has previously been presented for the phase- 
coherent correction of interstellar dispersion smearing. 
Convolution is performed by multiplying the spectrum of 
baseband data with the inverted frequency response of the 
interstellar medium (ISM), as modeled by the cold plasma 
dispersion relation (Hankins & Rickett 1975). The cur- 
rent development extends the concept of convolution to 
multiplication of the vector spectrum by the inverse of 
the frequency response matrix. The term "matrix con- 
volution" is used to distinguish this operation from two- 
dimensional convolution, such as would be performed on 
an image. Whereas the frequency response matrix may 
be composed of various factors, such as those arising from 
the ISM and ionosphere, this paper will be restricted to 
considering only the instrumental response. 

Following a brief review of matrix convolution and its re- 
lationship to the transformation properties of radiation, a 
technique is proposed by which the instrumental frequency 
response matrix may be determined and used to calibrate 
phase-coherent baseband data. As a preliminary illus- 
tration, the method is applied to full-polarization obser- 
vations of the millisecond pulsar, PSR J0437— 4715, from 
two different receiver systems. The results agree well with 
previously published polarimetry for this pulsar (Navarro 
et al. 1997). 

2. MATRIX CONVOLUTION 

Consider a linear system with impulse response, j{t). 
Presented with an input signal, e{t), the output of this 
system is given by the convolution, e'{t) = j{t) * e{t). In 
the two-dimensional case, each output signal is given by a 
linear combination of the input signals, 

e[{t)=jii{t)*ei{t)+n2(t)*e2{t), (1) 

e2(i) = J2i(i) * ei{t)+j22{t) * e2{t). (2) 
Now let ei{t) and 62 (i) be the complex- valued analytic 
signals associated with two real time series, providing the 
instantaneous amplitudes and phases of two orthogonal 
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senses of polarization. By defining the analytic vector, 
e{t), with elements ei{t) and 62 (i), and the 2x2 impulse 
response matrix, j{t), with elements jmn{t), we may ex- 
press the propagation of a transverse electromagnetic wave 
by the matrix equation, 

e'{t)=m*e{t). (3) 

By the convolution theorem, equation (3) is equivalent to 

E'iiy) = J{iy)E{i^), (4) 

where J(z^) is the frequency response matrix with elements 
Jmn{t^), and E{iy) is the vector spectrum. In the case of 
monochromatic light, or under the assumption that J{v) is 
constant over all frequencies, matrix convolution reduces 
to simple matrix multiplication in the time domain, as tra- 
ditionally represented using the Jones matrix. However, 
because these conditions arc not physically realizable, the 
Jones matrix finds its most meaningful interpretation in 
the frequency domain. 

The average auto- and cross-power spectra are neatly 
summarized by the average power spectrum matrix, de- 
fined by the vector direct product, P(z^) = {E{v)^E'^ {v)), 
where E^ is the Hermitian transpose of E and the angular 
brackets denote time averaging. More explicitly: 

p^^^ _ ( {EMEl{y)) {EME*2{y)) \ 
- \ {E^{u)El{u)) {E2{v)E*^{u)) ) ■ 

Each component of the average power spectrum matrix, 
Pmn (j^) , is the Fourier transform pair of the average corre- 
lation function, Pmn{T) (Papoulis 1965). Therefore, P(z^) 
may be related to the commonly used coherency matrix, 

p={E{t)®E\t)) = -p{<:)) = — P{v)dv, (6) 

where V{) is the centre frequency and 2Ai^ is the band- 
width of the observation. The average power spectrum ma- 
trix may therefore be interpreted as the coherency spectral 
density matrix and, in the limit Ai^ ^ 0, p = P(i^o)/27r. 

Using equations (4) and (5) it is easily shown that a two- 
dimensional linear system transforms the average power 
spectrum as _ 

P'(^)=J(:.)P(:.)jt(^). (7) 

This matrix equation is called the congruence transforma- 
tion, and forms the basis on which the frequency response 
of a system will be related to the input (source) and out- 
put (measured) coherency spectrum. For brevity in the 
remainder of this paper, all symbolic values are assumed 
to be a function of frequency, v, unless explicitly stated 
otherwise. 

3. ISOMORPHIC REPRESENTATIONS 

The subject of radio astronomical polarimctry received 
its most rigorous and elegant treatment with the two sepa- 
rate developments of Britton (2000) and Hamaker (2000). 
Whereas Britton illuminates the isomorphism between the 
transformation properties of radiation and those of the 
Lorentz group, Hamaker illustrates and utilizes the sim- 
ilarity with the multiplicative quaternion group. The 
salient features of both formalisms provide a sound founda- 
tion for the current development, and merit a brief review. 

Any complex 2x2 matrix, J, may be expressed as 

J = Jol + J • <T, (8) 



where Jq and J = (Ji, J2, J3) are complex, I is the 2x2 
identity matrix, and cr is a 3-vector whose components 
are the Pauli spin matrices. The 4- vector, [Jo, J], may 
be treated using the same algebraic rules as those used for 
quaternions. In the remainder of this paper, the equivalent 
quaternion and 2x2 matrix forms (as related by eq. [8]) 
will be interchanged freely. 

If J is Hermitian, then the components of [Jq, J] are real. 
The average power spectrum and coherency matrices are 
Hermitian and, when decomposed in this manner, 2[S'o. S] 
may be interpreted as the mean Stokes parameters, where 
2So is the total intensity and 2S is the polarization vector. 
By writing P = [6*0, S], it is more easily seen that the inte- 
gration of equation (6) leads to bandwidth depolarization 
when the orientation of S varies with v. 

An arbitrary 2x2 matrix may also be represented by 
the polar decomposition, 

J = JB^(/3)Rft(.^), (9) 
where J = (det J)l/^ 

Bm(/3) = exp(<7'm/?) = [cosh/3, sinh/3 rh], (10) 

P-ft('/') = exp(i <T' «(?!)) = [cos 0, z sin (/> n] , (11) 

and rh and n are real-valued unit 3-vectors. The matrix, 
Rft ((/)), is unitary and, beginning with equation (7), it can 
be shown to preserve the degree of polarization and rotate 
S about the axis, n, by an angle 2(j). Similarly, the Hermi- 
tian matrix, Bjfi(/?), can be shown to perform a Lorentz 
boost on the 4- vector, [So, S], along the axis, rh, by a 
velocity parameter 2/3 (Britton 2000). 

As both Bto,(/?) and Rft((/>) are unimodular, the congru- 
ence transformation (eq. [7]) preserves the determinant up 
to a multiplicative constant,_ | Jp,_ie. detP' = jJj^detP. 
Britton notes that detP = 5q — |5p = Sf^^ is simply the 
Lorentz invariant, which he calls the polarimetric invari- 
ant interval. He also proposes that a mean pulsar profile 
formed using the invariant interval may be used for high- 
precision pulsar timing. Whereas the boost component of 
the system response distorts the total intensity profile in 
a time-dependent manner, the invariant interval remains 
stable, providing a superior basis for arrival time estimates 
and a robust alternative to polarimetric calibration. How- 
ever, the invariant interval is not preserved in the presence 
of bandwidth depolarization, further motivating the cur- 
rent development. 

4. DETERMINATION OF THE SYSTEM RESPONSE 

In principle, any method by which the frequency re- 
sponse matrix may be derived with sufficiently high res- 
olution could be adapted for use with the matrix con- 
volution approach. A number of different parameteriza- 
tions and techniques have previously been presented for 
the calibration of single-dish radio astronomical instru- 
mental polarization (Stinebring et al. 1984; Turlo et al. 
1985; Xilouris 1991; Britton 2000). The current treatment 
follows the parameterization chosen by Hamaker (2000), 
including polar decomposition of the system response fol- 
lowed by application of the congruence transformation. 
This approach offers a number of significant advantages. 
For instance, in the derivation of the solution, it is unnec- 
essary to make any small-angle approximations and, when 
compared to the manipulations of 4x4 Mueller matrices, 
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the equivalent quaternion and 2x2 Jones matrix repre- 
sentations greatly simplify the required algebra. As well, 
equations (7) and (9) are specific to no particular basis, 
permitting the application of the formalism to a variety 
of feed designs by proper choice of the basis matrices, cr. 
Perhaps most importantly, as shown by Hamaker and dis- 
cussed below, the polar decomposition enables the deter- 
mination of the boost component of the system response 
using only an observation of an unpolarized source. This 
simplification has beneficial impact on experiments where 
only the total intensity or fractional degree of polarization 
are of concern. 

Consider the reception of unpolarized radiation, which 
has an input average power spectrum matrix, P^^I To/2, 
where I is the identity matrix and Tq is the total intensity. 
Beginning with equations (7) and (9), the output average 
power spectrum can be trivially shown to be 

P'i = B^ipf \J\^To/2. (12) 

Notice that the phase of J is lost in the detection of the 
average power spectrum. For this reason, the current tech- 
nique is insensitive to absolute phase terms arising in the 
frequency response of the system, and may be used to 
determine only the relative phase differences between its 
components. The scalar factors, |Jp and Tq, are easily 
determined in a separate flux calibration procedure that 
will not be considered presently. Therefore, Brfi(/3) may 
be found by taking the positive Hcrmitian square root of 
the measured average power spectrum matrix of an unpo- 
larized source (see Appendix A). 

This simple result merits closer inspection. The system 
response may be corrected by inverting equation (4), that 
is, by calculating E = 3~^E' . Notice that, under the polar 
decomposition chosen in equation (9), the boost transfor- 
mation is the first to be inverted. That is, regardless of 
the unknown rotation, the boost solution may be used to 
completely invert the distortion of total intensity. So. This 
consequence of the polar decomposition has considerable 
impact in the field of high-precision pulsar timing, where 
the average total intensity profile is used to determine the 
pulse time of arrival. Polarimetric distortions to the total 
intensity can significantly alter the shape of this profile 
and systematically alter arrival time estimates, especially 
as a function of parallactic angle, or with fluctuations of 
ionospheric total electron content. With the boost com- 
ponent thus determined, the resulting distortions to total 
intensity can, in principle, be completely corrected. 

It remains to solve for the rotation component of J, 
which may be determined using observations of calibra- 
tors with known polarization. Given the known boost, 
Brh,(/3)j input polarization state, P„, and measured out- 
put state, P^, equation (7) is solved for the rotation, 
R-n,(0)- Considering the equivalent three-dimensional Eu- 
clidian rotation, Rn^^\2(f>), of the polarization vector, S, 
it is easily seen that, given a single pair of input. Si, and 
output, S±' , polarization vectors, the rotations that solve 
Si' = Rft(3)(20)Si, where \Si'\ = \Si\, form an infinite 
set of rotations with axis confined only to a plane. There- 
fore, an observation of a second, non-collinear calibrator 
source is required in order to uniquely determine the sys- 
tem response rotation. 

Many receivers are equipped with a linearly polarized 
noise diode that may be used to inject a calibrator signal 



into the feed horn. This noise diode may be switched us- 
ing a wide-band, amplitude-modulated square-wave. The 
"off" or "low" fraction of the wave consists of only the 
system plus sky temperature, which shall be assumed to 
be unpolarized. The "on" or "high" fraction of the wave 
contains additional linearly polarized radiation, described 
in quaternion form by 

Pij = [I,(cos2*,sin2*,0)] Co/2, (13) 

where Co is the flux and ^ is the position angle of the 
calibrator diode. A single linear noise diode provides only 
one known input calibrator state. Unless another calibra- 
tor is available, the technique must be modified to solve 
for a reduced representation of the system response. For 
the calibration described in this paper, Rfi((/)) was be de- 
composed into two rotations: Ri)(A\E'), allowing imperfect 
alignment of the noise diode; followed by allowing 
a differential path length between the two linear polariza- 
tions (see Appendix A). It is important to distinguish A^E" 
and <!>/ from the actual parameters of the receiver, as the 
polar decomposition does not model the order in which 
the transformations physically occur. 

4.1. Calibrator Observations 

When observing the artificial calibrator, it is impor- 
tant that the noise diode be switched on a time-scale 
much shorter than the interval over which the digitization 
thresholds are reset, in order that P^ may be differenti- 
ated from the baseline, P l ■ However, the calibrator period 
must also be long enough to provide distinct on-pulse and 
off-pulse time samples in the synthesized filterbank. Given 
the sampling interval, ts, of the baseband recorder, the cal- 
ibrator period, Tc, should be at least Tq > ribtgN, where 
N is the number of channels in the synthetic filterbank and 
Ub is the desired number of phase bins in the integrated 
calibrator profile. In order to achieve clear separation of 
on-pulse and off-pulse states, rib = 64 was chosen. As well, 
it was found that a filterbank with N = 2048 channels suf- 
ficiently resolved the features of the frequency response. 

For each hour-long observation of PSR J0437— 4715, the 
pulsed calibrator was recorded for 4.5 minutes, ensuring 
sufficient signal-to-noise in each channel of the synthetic 
filterbank. Data were reduced offline using psrdisp, a 
software package developed to process pulsar baseband 
data, described in the following section. After forming a 
synthetic filterbank, the Stokes parameters in each channel 
were detected and folded at the pulsed calibrator period. 
High and low state polarimetric passbands were formed 
from the average on- and off-pulse Stokes parameters of the 
calibrator profile in each channel. These passbands were 
median filtered to remove spurious radio-frequency inter- 
ference, and interpolated to the required frequency resolu- 
tion, as dictated by the parameters of the coherent dedis- 
persion kernel. From this high frequency-resolution polari- 
metric representation of the calibrator, the frequency re- 
sponse matrix was computed as described in Appendix A. 
Figure 1 plots representative examples of the determined 
frequency response parameters at 660 MHz and 1247 MHz. 
The inverse of the frequency response matrix was used in 
the convolution kernel when reducing the pulsar observa- 
tions, as described in the next section. 
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5. PULSAR OBSERVATIONS 

Dual-polarization baseband data from the nearby mil- 
lisecond pulsar, PSR J0437— 4715, were recorded during 
2001 June 26-28, using the 64 m dish at the Parkes Ob- 
servatory. As a test of phase-coherent calibration in the 
pathological case (refer to Figure 2), the 50 cm receiver 
was used to record 3.6 hours of data at 660 MHz. As 
well, a total of 7 hours was recorded at 1247MHz, us- 
ing the centre element of the Multibeam receiver. The 
radio and intermediate frequency signals were quadrature 
downconvcrtcd, band-limited to 20 MHz, two-bit sam- 
pled, and recorded using the Caltech-Parkcs-Swinburne 
Recorder (CPSR) (van Straten, Britton, & Bailes 2000). 
The digitized data were processed offline at Swinburne 
University's supercomputing facilities using psrdisp, a 
flexible software package that implements a number of 
baseband data reduction options. 

The four channels of two-bit quantized data were cor- 
rected using the dynamic level-setting technique (Jenet 
& Anderson 1998) and combined to form the digitized 
analytic vector, e{ti). This vector was convolved while 
synthesizing a 16-channel coherent filterbank (Jenet et al. 
1997), a process described in more detail in Appendix B. 
In each channel, the Stokes 4- vector, [Sa,S]{i'k,tn), was 
detected and integrated as a function of pulse phase, 
given by t„ modulo the predicted topocentric period. 
Pulse period and absolute phase were calculated using 
a polynomial generated by the tempo software package 
(http: / /pulsar, princeton.edu / tempo) . 

Each 1 GB segment of baseband data (representing ap- 
proximately 53.7 seconds) was processed in this manner, 
producing average pulse profiles with 4096 phase bins, or 
a resultant time resolution of approximately 1.4 /xs. Each 
archive was later corrected for parallactic angle rotation of 
the receiver feeds before further integrating in time to pro- 
duce hour- long average profiles. These were flux calibrated 
using observations of Hydra A, which was assumed to have 
a flux density of approximately 85 Jy at 660 MHz and 48 Jy 
at 1247 MHz. The flux calibrated archives were then in- 
tegrated in time and frequency, producing the mean pulse 
profiles presented in Figure 3. In these plots, the linearly 
polarized component, _L, is given by = + U"^ and 
the position angle, ip, where tani/; = U/Q^ has been plot- 
ted twice, at i]j ± 7r/2. Comparison with the uncalibrated 
profiles in Figure 4 indicates significant restoration of the 
polarized component at both frequencies. As well, the cal- 
ibrated average profiles agree satisfactorily with previously 
published polarization data for PSR J0437— 4715 (Navarro 
et al. 1997). 

6. DISCUSSION 

Prior to its reception, the radiation from a pulsar must 
propagate through the magnetized plasma of the ISM, re- 
sulting in Faraday rotation of the position angle, tp. The 
change in position angle is given by A^/; = RMX^, where 
RM is the rotation measure and A is the wavelength. 
Therefore, one motivation for a pulsar polarimetry experi- 
ment is the determination of the RM along the line of sight 
to the pulsar. Combined with a dispersion measure esti- 
mate, the RM can provide information about the galactic 
magnetic field strength and direction. When the RM is 
very small, as it is for PSR J0437-4715, Ail; is not de- 



tectable across a single 20 MHz bandpass. Therefore, two 
observations widely separated in frequency are required 
in order to obtain an estimate. However, as there exists 
an unknown feed offset between 50 cm and Multibeam re- 
ceivers that remains uncalibrated in the current treatment, 
no estimation of rotation measure is presently offered. 

For pulsars with known RM, matrix convolution enables 
phase-coherent Faraday rotation correction by addition of 
Alp to A^ in equation (Al). This approach would further 
treat the problem of bandwidth depolarization, especially 
for sources with large rotation measures observed at lower 
frequencies. 

As phase-coherent calibration corrects the undetected 
voltages, any measurements derived from these data will 
also be implicitly calibrated. This includes not only the 
coherency products but also other statistical values such as 
the higher order moments, modulation index, and polar- 
ization covariance. The method may therefore find appli- 
cation in experiments which aim to describe pulse shape 
fluctuations (Jenet, Anderson, & Prince 2001) or varia- 
tions of the galactic magnetic field along the line of sight 
(Melrose & Macquart 1998), for example. 

Although matrix convolution may in principle be used 
to completely correct baseband data, a full solution of the 
system response is lacking in the current treatment. Per- 
haps a technique may be developed which combines the 
use of pulsars (Stinebring et al. 1984; Xilouris 1991) with 
high-resolution determination of the frequency response 
matrix. This approach would include tracking a bright, 
strongly polarized pulsar over a wide range of parallac- 
tic angles, followed by fitting a model of the instrumental 
response to the measured Stokes parameters. As it is non- 
trivial to model the effects of phenomena that may change 
over the course of such a calibration, such as fluctuations in 
ionospheric total electron content, it may prove more fea- 
sible to install a second artificial noise diode in the receiver 
feed horn. Oriented with a position angle offset from the 
first by approximately 45°, it would provide a more readily 
available and reliable means of regularly solving for Rji((^). 
Alternatively, a single noise diode could be mechanically 
rotated with respect to the feeds, obviating the inevitable 
cross-talk that would arise between two diodes. Consider- 
ation of these options may benefit future feed horn design. 

The observations of PSR J0437-4715 presented in the 
previous section serve to illustrate the strength the matrix 
convolution technique, highlighting its ability to recover 
data which has been severely corrupted by the observa- 
tory instrumentation. Perhaps the most rigorous test of 
the methodology will be derived from the unsurpassed tim- 
ing accuracy of PSR J0437-4715 (van Straten et al. 2001). 
Calibration errors translate into systematic timing errors, 
and any deficiency in the characterization of the system 
response would manifest itself in the arrival time residuals 
of this remarkable pulsar. 

The Parkos Observatory is part of the Australia Telescope 
which is funded by the Commonwealth of Australia for opera- 
tion as a National Facility managed by CSIRO. This research 
was supported by the Commonwealth Scholarship and Fehow- 
ship Plan. I am grateful to Matthew Britton and Simon John- 
ston, with whom I have shared many stimulating discussions on 
radio polarimetry. Thanks also to Matthew Bailes and Stephen 
Ord for assistance with observing and helpful comments on the 
text. 
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APPENDIX 



SOLUTION OF SYSTEM RESPONSE 

The instrumental frequency response matrix used in calibrating the data presented in this paper was parameterized by 

J = B^(/?)R^($j)R«(A*). (Al) 
The boost component may be solved most easily using the quaternion form of equation (12), 

[L'o, L'] = exp(2 cr. m/J)] J|^To = [cosh(2/?), sinh(2/3)Tfi] | J^Tq, (A2) 
where [Zq,l'] are the measured off-pulse Stokes parameters, producing 



, and rh=^T-. (A3) 

The two rotations, R^(4>/) and R{)(A^'), may be determined by considering the equivalent three-dimensional Euclidian 
rotations, R,5^^^(2$/) and R^^^' (2A\1/), of the input polarization vector. Given that the noise diodes installed in the 
receivers at the Parkes radio-telescope have a position angle nearly equal to 45°, the input Stokes parameters are [1, H]Co, 
where H = (0, 1,0) (cf. eq. [13]). Therefore, and A^f may be found by solving 

fl-" = Rq(3)(2$j)Rc(3)(2A*)fl-, (A4) 

where H" is the normalized polarization vector after the observed Stokes parameters have been corrected for the boost. 

It is given by 

[Ho, H]" = B^(/3)tp^B^(/3)/(| JI^Co), (A5) 
where are the observed Stokes parameters. Prom equation (A4), 

^ i„„-i / -^3 \ „„j A,T, _ 1 J.„„-1 / ^1 



= - tan-i — f and A* = - tan'^ -— — i— — . (A6) 

' 2 \H- J 2 V^2/cos(2$,); ^ ' 

Notice that | Jp, Tg, and Co cancel out in each of equations (A3) and (A6), obviating the need to solve for these parameters 
at this stage. The vector components of the boost quaternion, (Bi, B2, -B3) = sinh/3m, as well as and A^*, are plotted 
as a function of frequency in Figure 1. 



PERFORMING MATRIX CONVOLUTION 

It is assumed that the reader has some familiarity with the more common, one-dimensional form of cyclical convolution 
as it is performed in the frequency domain using the Fast Fourier Transform (FFT) (Press et al. 1992, §13.1). In the 
two-dimensional vector case, there arc simply two unique processes sampled at the same time interval, say p{ti) and q{ti). 
A one-dimensional, A'^-point, forward FFT is performed separately on each of p and q, forming two spectra, P{yk) and 
Q{i^k), < k < N. Corresponding elements from each of the spectra are treated as the components of a column 2-vector, 
E{vk) = {P{vk), Q{vk)), and multiplied by the inverse of the frequency response matrix, forming E' {vk) = J~^{i'k)E{i'k) 
(cf. eq. [4]). The components of the result, P'{vk) and Q'{vk), are once again treated as unique spectra, and separately 
transformed back into the time domain using the one-dimensional backward FFT. 

In the case of the pulsar data presented in this paper, p{ti) and q{ti) are the signals from the two linear feeds in the 
receiver, and 3{vk) consists of the instrumental frequency response matrix, as determined in Appendix A, multiplied by 
the dispersion kernel, H{vk) (Hankins & Rickctt 1975). When synthesizing an Af-channel coherent filterbank, H{vk) is 
divided in frequency into M distinct dispersion kernels, each tuned to the centre frequency of the resulting filterbank 
channel. 
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Fig. 1. — Frequency response at 50cm (top) and 20cm (bottom). For each band, the upper panel plots the rotation angles, A'l' (dashed 
line) and (solid line), and the lower panel plots the boost components along q (solid line), u (dotted line), and v (dashed line). A signal 
path length mismatch between the two polarizations from the Multibeam (20 cm) receiver results in a linear differential phase ("!>/) gradient 
across the bandpass. 




Fig. 2. — Average passbands from the two linear polarizations at 660 MHz. The gain in Ey (dotted line) varies wildly with frequency, 
evidencing the less than ideal condition of the 50 cm receiver. 
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Fig. 3. — Average polarimetric pulse profiles of PSR J0437— 4715 at 660 MHz (top) and 1247 MHz (bottom). The bottom panel of each plot 
displays the total intensity, Stokes / (upper solid line), the linearly polarized component, L (lower solid line), and the circularly polarized 
component, Stokes V (dotted line). The position angle, Vi in the upper panel of each plot is shown without respect to any absolute frame of 
reference. 
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Fig. 4. — Uncalibrated average polarimetric pulse profiles of PSR J0437— 4715 at 660 MHz (top) and 1247MHz (bottom). Severe differential 
gain variation in the 50 cm bandpass (see Figure 2) has decimated the linearly polarized component at 660 MHz (top plot, lower solid line), 
whereas the differential phase gradient across the 20 cm bandpass (see Figure 1) has nearly eliminated the circularly polarized component at 
1247MHz (bottom plot, dotted line). 



